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ABSTRACT 


We  evaluate  the  Bell-lattice-model  for  water  by  exploring  the  parameter 
space  by  means  of  the  Cluster  Variation  method  of  Kihuchi.  We  use  ten  order 
parameters  to  describe  the  different  occupation  of  tetrahedrons  in  the 
lattice,  with  and  without  hydrogen  bonds.  The  phase  diagrams  are  represented 
in  the  density-temperature  plane  and  in  the  pressure  temperature  plane. 

An  explanation  is  given  as  to  why  there  are  first  order  terminal  points  in 
the  metastable  region.  We  discuss  the  190®- rule  and  display  the  ratios  of 

hydrogen  bonded  pairs  versus  non-hydrogen  bonded  pairs. 
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1.  Introduction 

le  study  the  Bell*)  potential  for  the  lattice  model  of  water  using  the 
complete  cluster  variation  method.  The  basic  lattice  is  a  bcc  lattice  and  the 
sites  are  occupied  by  oxygen  atoms.  The  hydrogen  atoms  are  oriented  along 
the  lines  connecting  the  neighboring  atoms.  Since  the  bond  angle  between 
two  hydrogen  atoms  is  approximately  104° ,  the  hydrogen  atoms  should  lie  along 
those  two  body  diagonals  that  make  the  largest  angles.  This  is  accomplished 
by  assigning  appropriate  weight  factors.  The  complete  set  of  tetrahedral 
configurations,  as  well  as  all  the  sub-clusters,  is  given  in  reference  2.  We 
will  refer  to  this  model  as  JWAT  for  "Junior"  Water. 

Rather  than  use  the  10  clusters,  with  the  constraint  that  the  total 
probability  should  be  one,  we  prefer  to  introduce  the  nine  independent 
correlation  functions:  i.e.  all  completely  occupied  clusters  and 
subclusters.  (Table  1)  All  other  configuration  probabilities  can  be  written 
as  a  linear  function  of  these  above  defined  correlation  functions. 


There  are  18  configurational  probabilities.  The  factors  C  and  are 
given  in  Table  2.  If  the  configuration  is  empty  CA  iS  one,  otherwise  it  is 
zero,  the  coefficient  R  can  be  obtained  from  the  CR-table  in  reference  2,  by 
an  inversion. 

There  is  an  important  technical  point  to  mention.  In  the 
cluster-variation  method  one  can  express  each  subcluster  as  a  linear 
combination  of  larger  (sub)clusters.  For  instance,  if  on  a  given  point  the 
molecule  is  oriented  in  one  of  the  12  possible  ways,  this  will  determine 
whether  this  same  molecule  as  partner  in  a  pair,  can,  or  cannot,  be  hydroger 


bonded.  That  means  a  specific  choice  is  made  as  to  whether  the  linear 
combination  of  the  y's  will  contain  y4  (the  pair  without  H-bonding)  or  y5  (the 
pair  with  H-bonding).  In  the  JWAT  model  this  is  not  done,  instead  a  linear 
combination,  or  average,  of  these  two  options  is  taken.  This  average  is  taken 
with  certain  combinatorial  factors  in  accordance  with  the  so-called  ice  rules. 
*the  fact  that  we  did  not  specify  the  orientation  of  the  molecules  will  keep 
the  number  of  cluster  variables  low,  but  it  will  oversimplify  the  entropy 
expression  and  it  will  of  course  bar  us  from  any  information  about  the 
possible  rotational  ordering  of  the  molecules. 

A  second  simplification  we  use  in  the  JWAT  -  model  is  the  half-filled 
state,  i.e.  density  approximately  one  half.  In  this  model  we  find  this  state 
at  low  temperatures  by  adjusting  the  energy  parameters.  Such  a  state  refers 
to  preferred  occupation  of  those  clusters  that,  for  a  given  range  of  values 
for  the  chemical  potential,  have  two  sites  occupied.  It  is  not  specified  as 
to  which  two  out  of  four  sites  are  occupied,  the  model  again  takes  an  average 
over  two  options.  2.  Interaction  parameters 

The  lattice  model  contains  three  parameters: 

the  pair  interaction  between  nearest  neighbors  when  they  are 
not  hydrogen  bonded. 

eH  the  additional  pair  energy  between  two  neighboring 

molecules  when  they  are  hydrogen  bonded. 

ennn  repulsive  next-nearest  neighbor  interaction. 

The  first  two  energies  are  attractive;  the  symbol  stands  for  the 
absolute  value  of  a  negative  quantity.  All  energies  are  expressed  in  units  kR 
and  hence  have  the  dimension  of  a  temperature. 
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The  absolute  values  of  the  interactions  are  only  important  for  scaling 
purposes.  Using  the  principle  of  corresponding  states  the  topology  of  the 
phase  diagram  depends  only  on  a  reduced  set  of  interaction  parameters. 
Generally  one  divides  by  the  largest  (i.e.  mostly  the  nearest  neighbor 
interaction)  parameter.  In  the  case  of  the  water-model  one  has  the  choice 

between  ep>  or  ^  (35  Dell  used)  or  ep  +  e^,  the  total  binding  energy  of  a 
hydrogen  bonded  pair.  We  prefer  the  first  choice. 

In  fig.  1  and  2  the  potential  diagrams  for  the  first  two  cases  are 
drawn,  together  with  the  values  and  mnemonics  of  the  actual  potentials  for 
which  computations  have  been  performed.  The  precise  values  are  collected  in 
"Bible  3.  Potential  10  is  the  same  as  potential  16  without  hydrogen  bond 
energy.  Potentials  16  and  17  have  been  treated  in  an  extended  mean  field 
model  by  Bell  and  Shit.3)  (our  potential  16  is  slightly  different  due  to 
rounding  errors;  Bell  and  Salt  have  ©p/e^  =  1.45  and  we  have  ~ 

1.450666...) 

There  is  a  general  limitation  on  the  choice  of  potentials,  determined  by 
our  request  that  the  intermediate  state  exists.  We  calculate  the  generalized 
pressure  from  ft  ^ 

and  for  T  ■  0  there  are  three  states  possible:  ^  =  0  (vacuum)  with  P  =  0,  (3 

“  1/2  (open  ice)  with  P  +  1/2  p  +  ep  +  and  ^sl  (closed  packed)  with  P  +  JJ 

+  4ep  +  2  eh  -  3  If  we  plot  P  at  T*  0  (compare  fig.  3)  find 

three  possible  transitions  at 

p0  —  —2-  t  l 

IPt  =  -V  to  -w.) 
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In  order  to  obtain  a  stable  intermediate  phase  we  need 


=~  -€f>  +  %eM»u 

this  leads  to  the  "forbidden"  area  in  figures  1  and  2.  The  values  of 

are  given  in  Thble  3. 


>s° 

P«, 


The  value  of  the  hydrogen  bonding  parameter  does  not  come  in  the  above 
condition,  since  we  compare  perfect  open  ice  with  the  close  packed  phase.  The 
energy  difference  between  c.p.  and  open  is  given  by  two  more  nearest  neighbor 
bonds  and  three  more  next-nearest  neighbor  bonds  as  can  be  found  from  an 
inspection  of  the  lattice  structure  when  it  is  either  half  occupied  (diamond 
lattice)  or  fully  occupied  (bcc).  The  hydrogen  bond  energy  does  come  into 
play  in  the  P(p)  lines  for  each  of  the  regions  below  and  above  ,  i.e. 
is  in  the  open  and  in  the  c.p.  section  separately. 


3  The  entropy  at  zero  temperature 


Since  in  the  JWAT  model  no  sublattices  are  introduced,  and  no 
orientational  effects  are  considered,  except  for  the  hydrogen  bonding,  the 
open  and  closepacked  ices  cannot  be  treated  completely  correctly  by  this 
model.  In  order  to  estimate  the  difference  in  entropy  at  low  temperature 
between  the  various  models  we  computed  the  entropy  for  c.p.  and  open  in  JWAT 
.  The  open  ice  structure  in  the  JWAT  model  can  be  associated  with  the  case 
that  all  tetrahedra  are  in  the  configuration  4  in  reference  2  (or  24  in  Thble 


This  cluster  has  multiplicity  72:  4  due  to  the  distributions  of  the  hydrogen 

bond  on  the  four  nearest  neighbor  edges  of  the  tetrahedron,  and  I  8  because 


ftige  7 


of  the  orientational  multiplicity  of  the  hydrogen  bond. 

The  correlation  functions  of  the  subclusters  are  now: 

Pl6  -  1/48;  p20  =  1/72;  p12  =  1/24;  pn  =  1/4;  p12  = 

1/48;  P3  =  1/72;  p10  =  1/2  and  pt  =  1/24. 

Substituting  this  in  the  expression  for  the  entropy  in  units  k^; 

-  S  =  2.  30 

where  ^  are  the  Kikuchi  numbers,  ^  \  the  (total)  weight  factors  and  the 
sum  is  over  all  contributing  clusters.  The  result  is: 

= l>h]  3  a  ~  -2-7^2- 

The  close  packed  state  results  if  all  tetrahedra  are  in  configuration  8  in  ref 


2  (or  number  9  in  table  1). 


the  value  of  the  subclusters  is  now:  jp^  -  \  =  W  llib> 

’/^and  Pj=l/2.  This  leads  to  an  entropy 

=  -la  ^2.  4  ^7^ 

The  absence  of  a  proper  two-sublattice  description  in  the  JWAT-model,  we 
obtain  a  rather  large  negative  value  of  the  entropy  of  these  two  phases.  The 
result  is  that  the  pressure,  as  given  by 

—  p  ci  -  S 


will  also  be  much  lower  than  it  would  be  if  the  entropy  was  given  correctly. 
Hence  the  open  and  close  packed  ice  phases  are  underestimated  in  this  model  , 
compared  to  the  liquid/vapour  -  i.e.  the  disordered  phases.  This  is  most 
clearly  seen  by  comparing  the  two  sublattice  phase  diagrams  for  potentials  12 


and  17  with  their  JWAT  counterparts 


It  is  furthermore  interesting  to  compare  the  ground  state  entropy  with 
the  results  obtained  by  Bell.  Starting  with  his  equation  3.4,  i.e. 

le 

5?  -  -  r  -  t  g  i  Vc’1 

where  ‘UK  are  representing  the  ten  tetrahedrons,  their  weight  factors  and 
eo  given  by: 

^  c,  ■=  S  NJ  ^  9  +  (■ 1  1  ~  ?  H 


Translated  in  our  language,  he  took  for  the  Kikuchi  numbers  the  values:  1  for 
the  tetrahedrons  and  -3  for  the  points,  (and  all  others  zero).  That  is  he 
considered  a  collection  of  N  te trailed rons  which  overlap  only  with  their 
vertices,  hence  /f  NT  points  are  substracted  from  the  N  point  terms.  The 
systematic  procedure  of  the  CVM  is  to  consider  all  N  tetrahedrons  that  fill 
the  space  completely.  These  overlap  with  triangles,  second  neighbors,  first 
neighbors  and  point  sites  and  thus  is  taken  into  account  by: 

S  =  t  -IZ  T  -  % 


rather  than 


S  -  Sfr  -a  £ 


as  in  Bell.  It  turns  out  that  the  large  negative  contribution  comes  from  the 
triangular  term  -12  S^,.  For  comparison  we  give  also  the  numerical  values  of 


the  entropy  as  obtained  by  Bell's  expression. 


5u.  =  +  ^  -  f.s"5 

St-U  =■  -'i-Cx-'Z.-f  ^  ~  0.<j3 
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Finally  we  would  like  to  mention  that  the  Pauling  estimate  for  open  ice  leads 
to  So  =  1/2  In  3/2  =  0.2027  and  twice  that  amount  for  the  c.p.  state. 

4.  Description  and  Comparison  of  a  Number  of  Potentials;  Densi ty- 
Temperature  Phase  Diagram. 

If  we  look  at  the  0-T  plot,  the  influence  of  the  intermediate  state  will 
show  up  in  two  different  ways.  Either  we  find,  at  high  temperatures  a  direct 
transition  between  the  disordered  low  density  state  to  the  close  packed  state, 
or  we  find  a  transition  from  the  low  density  state  to  the  intermediate  state 
and  subsequently,  at  somewhat  lower  temperatures,  a  transition  from  the 
intermediate  density  state  to  the  close  packed  state.  To  make  an  easy 
distinction  between  these  two  case  we  will  refer  to  the  first  case  as  the 
one-dcme  phase  diagram,  i.e.  ^)-T  diagram,  and  to  the  second  as  the  two-dame 
diagram. 

A  second  characteristic  of  the  ^-T  plots  is  the  appearance  of  a  "slit",  a 
region  near  one  half  density  where  the  <^>  is  zero  to  one  transition  is 
interrupted  by  an  intermediate  state  so  that  the  transition  now  takes  place  in 
two  steps,  one  from  density  near  zero  to  a  density  somewhat  less  than  a  half 
and  one  from  a  density  somewhat  more  than  a  half  to  a  density  almost  one. 
Although  this  two  step  process  is  similar  to  the  two  step  process  describod 
above  in  the  "two  dome"  diagram,  there  is  an  important  distinction,  i.e.  at 
high  temperatures  this  transition  disappears.  The  point  of  disappearance  is  a 
triple  point  since  three  phases  are  in  equilibrium  with  each  other.  We  will 
see  later  that  for  certain  values  two  triple  points,  since  the  above  mentioned 
"slit"  does  also  occur  in  certain  two-dome  p-T  diagrams. 
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4,1)  Ttoo-dcme  potentials. 

Let  us  start  with  potentials  12,  15  and  11  that  have  a  relatively  large 

ennn/ep  a  variable  amount  of  ej{.  In  potential  12  we  have  a  two-dome  (See 
Pig.  4)  structure,  i.e.  two  critical  densitites.  The  low  density 
critical  temperature  is  at  400  K,  the  high  density  critical  temperature 
at  350  K.  The  lower  value  of  the  second  critical  temperature  is  due  to 
the  fact  that  at  these  density  values  the  fraction  of  nearest  neighbor  bonds 
that  is  hydrogen  bonded  is  less,  leading  to  a  analler  effective  pair 
interaction. 

Tb  obtain  the  high  density  transition  one  needs  a  higher  presure  and 
compared  to  low  density  critical  temperature  this  is  clarly  seen  in  the  P-T 
plot  for  this  potential  (Fig.  5).  from  this  plot  it  is  also  clear  again  that 
the  intermediate  phase  is  open-ended.  The  upper  branch  descends  and 
subsequently  rises  somewhat  before  it  ends  in  the  second  critical  temperature. 
This  second  terminal  point  causes  us  to  reject  this  type  of  phase  diagram, 
since  it  does  not  reconcile  between  this  type  of  result  and  even  a 
qualitative  comparison  between  this  type  of  result  and  the  observed  diagram 
for  water. 

A  similar  set  of  figures  (Fig.  6  and  7)  was  obtained  for  potential  15. 
At  low  temperatures  it  takes  considerably  more  pressure  to  accomplish  the 
second  transistion. 

The  third  potential  we  want  to  discuss  in  this  section  is  number  11. 
This  is  still  a  potential  with  two  domes,  but  it  has  a  characteristic  we  will 
find  in  many  of  the  one  dome  potentials:  there  is  a  triple  point  at  Tg. 
Below  Tg  the  potential  develops  a  "slit",  as  described  above.  Tg  lies  at 
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152.8  K.  The  composite  sketch  (Fig.  8)  of  the  T-Oresults,  gives  a  second 

|  ^ 

triple  point  at  lower  temperatures.  Below  this  temperature  we  have  two 
transitions  again.  In  the  intermediate  region  T^<T<T3  We  have  a  direct 
transition  from  the  low  to  the  high  density  state.  All  dashed  lines  are 
metastable  transitions.  Please  note  that  the  low-to-high  transition  extends 

above  T3>  but  ends  in  a  metastable  end  point  TM.  The  accompanying  pressure 
temperature  plot  is  given  in  Fig.  9.  The  lower  branch  1  is  related  to  T*, 
the  upper  branch  2  is  related  to  The  metastable  branch  that  goes  up  to  TM 
is  just  visible  above  branch  number  one.  The  two  triple  points  are  on  the 
left. 


Finally,  we  describe  potential  18,  which  is  almost  a  one-dome  potential. 
The  density  temperature  diagram  is  given  in  figure  10A  qualitatively  and  in 
figure  10B  quantitatively,  but  incomplete.  The  noteworthy  difference  between 
figures  10  and  8  is  that  the  metastable  states  between  the  two  triple  points 
are  no  longer  connected. This  means  in  the  P-T  plot  (see  figure  11),  that  the 
second  branch  now  consists  of  two  independent  parts,  one  labelled  "open/cp" 
and  another  labelled  liquid  (1/2)  -  liquid  (1).  This  last  line  extends  a 
little  bit  below  the  "vap./liq."  line. 

At  the  end  of  this  sub-section  we  would  like  to  discuss  the  termination 
of  phase  lines  in  the  P-T  diagrams.  As  a  rule  phase  lines  can  only  end  in 
second  order  critical  points  (ref.  4)  since  the  double  tangent  construction 
at  temperatures  just  below  will  relate  two  values  of  the  order  paramenter 

say  ji  and  ^  we  approach  Tc , ^  moves  towards  ^ 3  and  the  two  will 
coincide.  Not  so  for  a  metastable  double  transition.  Imagine  a  function  of 
at  least  the  sixth  order,  chosen  synmetric  for  convenience.  In  fig.  12a  we 
took  the  central  value  larger  than  the  two  maxima  at  the  left  and  the  right. 


We  can  make  two  double  tangent  construction  ,  leading  to  a  transition  from 

^ l>  to  ^2  a  second  transition  from  ^  to^  This  is  the  situation  below 
the  triple  point.  At  the  triple  point,  the  three  maxima  are  of  equal  height 
and  the  tangents  coincide:^  =  ^3.  Above  (below)  the  triple  point  we  deal 
either  with  a  direct  transition  from  ^  to  £>4  or  with  two  metastable 

transitions  from  ^  to  ^2  and  from  ^to  (^4,  but  ^  2  and  ^  3  appear  in  the 
reversed  order,  i.e.:  the  tangents  cross.  This  option  eventually  may 
disappear  (see  fig.  12d)  when  the  second  derivative  in  ^  ancj  ^  becomes 
zero.  This  is  the  condition  for  the  terminus  of  the  first  order  metastable 
transition.  Hence  it  is  incorrect  to  conclude  that  if  the  P-T  line  ends  one 
is  necessarily  dealing  with  a  second  order  transition,  at  least  not  for  a 
meta stable  branch. 

4.2  One  dome  potentials 

Potentials  16,  17,  13  and  14  are  all  of  the  "one-dome"  type,  that  is 
there  is  only  one  critical  point  at  high  temperatures  and  we  find  a  triple 
point  at  lower  temperatures,  which  leads  to  two  phase  transitions  at  a  given 
temperature  below  this  triple  point.  In  the  P-T  diagram  we  have  a  stable 
vapor-open  ice  coexisting  line  and  a  stable  open  to  close-packed  line.  Open 
ice,  the  state  where  the  density  is  approximately  one  half,  disappears  above 
the  triple  point  and  the  only  stable  coexistence  line  which  is  left  is  the 
vapor  to  close  packed  line,  which  in  this  simple  model,  is  indistinguishable 
from  the  vapor-liquid  line. 

Empirically  we  found  that  the  range  of  values  of  ennn  that  leads  to  this 
topology  is  found  in  the  area  of  figure  1  below  a  line  given  by 

^n»«/Gb  ~  (Cw/er) 
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and  above  the  horizontal  line  that  guarantees  three  groundstates,  as  discussed 
In  section  2. 

The  general  features  of  the  one-done  potential  can  be  seen  from  16  and 
17,  which  are  almost  coinciding.  See  fig.  13  for  the  temperature-density 
plot.  This  plot  does  not  show  that  there  is  a  triple  point  at  about  30. 
Figure  14  is  the  P-T  plot.  The  same  remark  holds  here;  there  is  a  branch  at 
temperatures  below  150,  the  left  axis  of  the  figure. 

Ratential  13  is  similar  to  16  and  17  but  since  e^  iS  larger  the 
rectilinear  diameter  is  more  bend  to  the  left  as  can  be  seen  in  figure  15.  It 
also  has  a  triple  point  but  at  lower  temperatures.  In  potential  14  the  triple 
point  lies  higher  with  respect  to  the  critical  point.  In  this  potential  we 
haVe:  Tc  =  480  311(1  T3  =  240«  866  fig111-®  The  relative  proximity  of  Tc 
and  t3  makes  it  inviting  to  compare  this  potential  with  real  water.  The 
pressure-  temperature  plot  (figure  17)  shows  an  undercooling  branch  far  into 
the  negative  pressures.  It  is  strongly  indicated  from  preliminary  reports 
about  the  work  of  R.  Speedy  5  gaich  a  curve  indeed  exists  in  water  that 
is  under  tension. 

4.3  Evolution  of  one  type  into  the  other 

Let  us  focus  on  potential  18,  figure  10A.  This  diagram  is  the  lynchpin 
of  our  considerations.  By  imagining  changes  one  can  obtain  all  others.  At 
one  hand;  if  the  two  triple  points  move  towards  each  other,  first  we  obtain 
figure  8  (no  metastable  terminal  points  in  the  center  section)  and  eventually 
figure  4  a  pure  two-dcme  phase  diagram,  because  the  two  triple  points  will 
eliminate  each  other,  as  can  be  seen  from  a  succession  of  P-T  curves.  At  the 
other  hand  if  the  upper  triple  points  moves  higher,  the  two-dcmc  feature 
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disappears  and  we  obtain  a  single  critical  temperature. 

5.  Behavior  around  a  triple  point.  A  triple  point  in  the  P-T  plot  means 

that  three  branches  meet  in  one  point.  Such  encounters  are  subject  to  "The 
180°-Rule"  (see  ref.  6).  Each  line  has  a  metastable  extension  given  by 

a  dotted  line  in  fig.  18  and  one  way  to  express  the  rule  is  that  the 
extension  of  line  (should  always  lie  between  the  lines  2  and  3.  Qir 
computations  satisfy  this  requirement. 

6.  The  hydrogen  bonding 

The  subcluster  variables  y^  to  y^  denote  the  fraction  of  empty  nearest 
neighbor  pairs,  the  fraction  of  nn  pairs  that  are  half  occupied  and  half 
unoccupied,  the  fraction  of  pair  sites  that  are  both  occupied  without, 
respectively  with,  hydrogen  bondingg.  It  makes  more  sense  t  o  consider  either 
the  ratio  between  hydrogen  bonded  pairs  and  bonded  pairs:  y3/x3. 

Plotting  of  the  first  ratio  is  of  some  advantage  since  it  can  easily  be 
checked  for  its  assumption  value  at  low  densities: 

JU. _  y«,„  _  ceH/k.rj__ 

7+  fGw/u 

since  the  ration  y^/y3  approaches  the  appropriate  Boltzmann  factor.  At  higher 
densities  this  approximation  f  ails  and  as  a  result  of  the  ice  rules  the  number 

of  hydrogen  bonded  pairs  y^  =  compared  to  the  number  of  non-hydrogen  bonded 
Ym  =  y3  starts  to  decrease.  Ekch  separately  does  increase,  at  higher 
densities,  but  y^  increase  more  rapidly. 

At  higher  densities  and  lower  temperatures  the  oxygens  will  make  as  many 
hydrogen  bonds  as  possible,  at  the  cost  of  the  non-hydrogen  bonded  pairs,  but 
eventually  this  can  no  longer  be  done  in  the  dissordered  phase  and  this  leads 
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to  instability. 


Numerical  results  for  the  ration  (here  plotted  as  yN/y{J)  are  given  for  a 
"two  dome"  potential  12  (see  figure  19)  and  for  a  "one  dome"  potential,  nunber 
14  (see  figure  20).  The  last  case  illustrates  the  description  we  just  gave: 

The  number  of  y^  at  very  high  densities  decreases,  but  at  low  temperatures 
this  phase  is  no  longer  stable.  The  stable  part  is  near  the  half  density,  in 
the  center  of  the  figure. 


We  also  recall  the  remark  made  in  the  first  paragraph  of  section  4.1, 
i.e.  densities  around  3/4  the  number  of  H-bonds  is  less,  than  at  de^nities 
near  1/4,  due  to  poor  "accomodation"  of  the  orientation,  which  leads  to  a 
lower  critical  temperature  at  the  higher  density.  The  most  sensible  plot  is 
yH/5  ,  but  the ^6  are  not  yet  available. 

7.  Phase-diagram  at  low  temperatures. 


In  this  section  we  consider  the  behavior  of  the  phase  diagram  at  low,  but 
non-zero,  temperatures,  in  particular  potential  14  (see  figures  16  and  17). 
There  is  an  almost  horizontal  line  in  the  P-T  diagram.  Below  this  line  we 
have  the  vapor  phase,  above  the  line  we  have  either  the  open  phase  (at  low 
temperatures)  or  the  close  packed  phase  (at  temperatures  above  approximately 
1X240).  The  dividing  line  starts  at  T=0  and  P=260,  the  "equality"  pressure 
mentioned  in  section  2.  This  line,  the  cp-open  transition  descends  towards 
the  horizontal  vapour  line  and  ends  at  a  triple  point  Pg,  Tg.  The  extension 
below  the  axis  is  a  metastable  transition. 


The  descending  P-T  curve  is  determined  by: 


dp/  _mp|  i-^=q 


Using  the  values  for  the  entropy  computed  in  section  3  gives 

n  «  /"fc>P  \  _  n  * 


This  is  not  the  case  in  the  two  orientation-two  sublattice  model,  we  will 
discuss  in  the  next  report.  The  negative  slope  is  due  to  the  fact  that  S  open 
<  1/2S  . 

In  the  origin  we  have  actually  two  lines  starting  at  P=0:  Che  is  the 
vapor  close  packed  which  is  metastable  up  to  the  triple  point  and  one  that  is 
related  to  the  vapor-open  transition.  The  last  ends  at  the  triple  point,  at 
least  its  stable  part,  and  carries  towards  that  point  with  a  different  slope 
(180°-rule).  A  simple  calculation  confirms  this.  For  vapor-open  we  have 


^  =  -."T  Cl ^  §  l 

At  T=0:  P  vap-open  “  Q  and  ^  vap_p  =  pi*  The  inequality 

p0  -L.  pj  leads  to 

Hp  y  ,  .  .  C 


'Vj  -  O^JL^ 

they  become  equal  at  the  triple  point. 
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Tables 

1*  The  nine  independent  correlation  functions  of  the  model,  expressed  as 
clusters  and  subclusters. 

2.  The  remaining  clusters  and  subclusters  expressed  as  linear  combinations  of 
the  nine  independent  correlation  functions. 


3.  Potentials  used  in  this  report,  unreduced  and  reduced  values  of  the 
parameters. 


Ffege  18 


Figures 

1  The  potentials  analyzed  in  this  report:  reduced  by  means  of  Ep. 

2  Same,  reduced  by  E^ 

3  Pressure  as  function  of  the  chemical  potential  at  T=0. 

4  Density-temperature  plot  for  potential  12. 

5  Pressure- temperature  plot  for  potential  13 

6  Density-temperature  plot  for  potential  15. 

7  Pressure-temperature  plot  for  potential  15. 

8  Density-temperature  plot  for  potential  11. 

9  Pressure- temperature  plot  for  potential  11. 

10  Density-temperature  plot  for  potential  18. 

11  P-T  plot  for  potential  18. 

12  Double  tangent  construction  and  metastable  states. 

13  Density-temperature  plot  for  potential  16  and  17. 

14  Pressure-temperature  plot  for  potentials  16  and  17. 

15  Density-temperature  plot  with  rectilinear  diameter  as  dashed  line 
for  potential  13. 

16  Density-temperature  plot  with  rectilinear  diameter  for  potential 
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17  Pressure-temperature  plot  for  potential  14. 

18  Illustration  of  the  180°-rule. 

19  The  ration  of  normal  bonds  to  hydrogen  bonds  (Potential  12). 

20  The  ratio  of  normal  bonds  to  hydrogen  bonds  (Potential  14). 
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